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As is well known, primordial gravitational waves may be amplified to detectable levels by parametric 
amplification during eras when their wavelengths are pushed outside the cosmological horizon; this can 
occur in both inflationary and "pre-big-bang" or "bounce" cosmologies. The spectrum of a gravita- 
tional wave background is expressed as a normalized spectral energy density Q(uj) = (oj / p c )(dp sv / did) , 
where p c is the critical energy density that makes the universe spatially flat, and dp gw is the energy 
density of gravitational waves in a frequency band duo. The logarithmic slope of Q is simply related 
to three properties of the early universe: (i) the gravitons' mean initial quantum occupation number 
N(n) (= 1/2 for a vacuum state), where n = au) is the (invariant) conformal frequency of the mode 
and a is the cosmological scale factor, and (ii) & (iii) the parameter 7 = p/p of the cosmological 
equation of state during the epoch when the waves left the horizon (7 = 7 i) and when they reentered 
(7 = 7/). In the case of an inflationary cosmology, the spectral index is 
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and for bounce cosmologies it is 
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These expressions are compared against various more model-specific results given in the literature. 



I. INTRODUCTION 

The prospect of a detectable cosmological background 
of gravitational waves has opened up a new avenue for 
the investigation of fundamental physics and cosmology. 
Gravitational waves couple very weakly to matter, and 
are not blocked or thermalized during any cosmological 
epoch back to Planckian densities. The waves do, how- 
ever, couple to large-scale cosmological spacetime curva- 
ture, and can therefore provide information about the 
evolution of this curvature, and hence about the funda- 
mental physics of the matter fields driving this evolution. 

A great deal of work has been done on this subject, 
deriving gravitational-wave spectra from various cosmo- 
logical models, usually in hope of producing a model 
that would generate a detectable signature in planned 
gravitational- wave detectors. The most comprehensive of 
these is a recent paper by Gasperini JjJ] , which presents a 
general prescription for constructing gravitational-wave 
spectra from fairly arbitrary cosmologies. The results in 
the present paper are consistent with Gasperini's results, 
although I derive them using a different formalism and 
consider a more generic cosmology. 

The primary intent of this paper, however, is more the 
reverse of this: rather than deriving a spectrum from 
a particular cosmological model, this paper shows how 
generic properties of the early universe can be imme- 
diately deduced from an observation of a gravitational- 



wave background in any spectral band. To this end, I 
have chosen as my observable the spectral characteristic 
that is least susceptible to the tunings of particular cos- 
mological models — namely, the spectral index — and 
have expressed it in terms of quantities that are most di- 
rectly related to the underlying physics of the cosmologi- 
cal "fluid" — namely, the equations of state p = 7 p. Fur- 
thermore, since the initial state of the gravitational waves 
(before amplification) is potentially one of the more in- 
teresting characteristics that might be deduced from ob- 
servations, I have left it as a free parameter, rather than 
making the usual assumption of starting the waves in a 
quantum ground state. 



A. Organization of this paper 

In Sec. ||, I present a skeleton of the derivation of the 
main formulae of this paper. Much of the derivation has 
been done in one form or another in the published lit- 
erature, and those familiar with the field will find noth- 
ing surprising, except perhaps the consideration of non- 
vacuum initial conditions in Sec. 



[IB. Since the result 



is largely a generalization of the more model-specific for- 
mulae found in the literature, Sec. Ill compares this 
per's formulae with those published previously. Sec 
presents some concluding remarks. 
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II. THE COSMOLOGICAL AND WAVE 
EQUATIONS 

I consider gravitational waves that are linear pertur- 
bations on the spatially flat FRW metric: 

ds 2 = a 2 (ri)[-dr] 2 + dx 2 + dy 2 + dz 2 ] + h a pdx a dx fi , 

(2.1) 

where the conformal time coordinate 77 is related to 
proper time t by dt = a{rf)drj (I have chosen units in 
which c = 1). The cosmology has a perfect fluid source 
that obeys the instantaneous equation of state p = jp, 
where 7 need not be constant but typically evolves slowly 
(7V7 *C a' /a, where ' = d/drj). Causality considerations 
require that 7 < 1; realistically, we can assume 7 < 1, 
since not all of the energy in the universe will be in the 
maximally-stiff field. One can solve Einstein's equations 
to zeroth order in h to obtain the exact solution 



a{r}) = do exp 



dm 



a /a' Q + j; i y m [l + 37(7? 2 )]/2 j ' 



(2.2) 



which reduces to the usual power-law evolution 0(77) ~ 
772/(1+37) during epochs of (nearly) constant 7. 

Following the notation in Eq. (2) of 0, I write the 
linear perturbations in the form: 



h a - V ^-U (x)e U) 
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where e^l is some basis of polarization tensors, J7 n (x) oc 
e m x is a spatial harmonic function with a true (physical) 
wave number k = n/a and frequency uj = n/a (n = ||n||), 
and /J, n (v) obeys the Schrodinger-likc equation: 



Mr, 



Mn = 
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Physically, the relative magnitudes of n 2 and the effective 
potential a" /a are related to whether the wave is inside 
or outside of the Hubble radius. Comparing k to the 
Hubble radius rji = a/(da/dt), one has: 
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Roughly speaking, a wave that is hitting the effective 
potential (a" /a increasing to meet n 2 ) is one that is exit- 
ing the Hubble radius, to order of magnitude. Similarly, 
a wave that is emerging from the effective potiential is 
reentering the Hubble radius. The exception is for ep ochs 
when 7 — ► 1/3; see the remarks near the end of Sec. [I A. 

Cosmological amplification of waves occurs where ini- 
tially oscillatory waves hit the effective potential during 



epochs of accelerated collapse (a' < 0, 7 > —1/3) or ac- 
celerated expansion (a' > 0, 7 < —1/3), when \a"/a\ is 
increasing, and then emerge into a post-inflationary uni- 
verse (a' > 0, 7 > —1/3) in which \a"/a\ is decreasing. 
A schematic of such an evolution is shown in Fig. [j]. The 
periods when the wave initially hits the potential and 
when it finally emerges from it will be denoted by sub- 
script i and /, respectively. The (complex) amplification 
factor of the emerging waves /3 n = u n (rj f) / u n (r]i) is 
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(2.6) 



this is Eq. (11) of [|| with a sign error corrected, and 
reexpressed in the current notation. 
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1. Schematic of the evolution of the effective po- 



FIG. 



tential \a" /a\ (sollid line) and Hubble scale (a'/a) 2 (dotted 
line) for inflation (77 < r/i) followed by decelerating expansion 
(77 > t?i ). Also shown are the times r]( and 77/ when a wave 
with squared conformal frequency n 2 (dashed line) hits and 
emerges from the potential barrier. 

If the evolution of a(rf) between at and a/ is monotonic, 
one can follow || and use the equation of state to reex- 
press the integral over 77 as an integral over a, eventually 
obtaining, for af 3> af. 
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(2.7) 



where (7) is the average of 7 from ?7i to 77/ with weighting 
factor 1/a 2 . 



The derivation of Eq. (2.7) breaks down for the case of 
a bounce cosmology , wh ere 0(77) is not monotonic. In this 
case, however, Eq. (2.6) is almost entirely dominated by 
the contribution of the integral J dr\j a 2 near the time of 
the bounce. Modeling the bounce as smooth with some 
finite concavity <Zq around the minimum eto, we obtain: 
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A. Dependence on n 

The amplification factor (3 n is at least implicitly de- 
pendent on n, since waves of different n will hit and 
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emerge from the effective potential at different times 
rji and rjf, and hence with different aj, aj. For peri- 
ods of constant 7, we have a ~ ^ 2 /( 1 + 3 7) > The con- 
dition n 2 = a" /a for hitting or emerging from the ef- 
fective potential implies that aij ~ n~ 2 /( 1+37i >^, and 

constant. If 7 is 



that a' i f /(a,ijn) — \/2/|l - 3^,/ 1 
changing gradually near j^y these formulae remain true 
up to corrections of order (j'/^)(o/Ja) . Furthermore, 
the weighted average (7) in Eq. (2/7) depends only very 
weakly on rjij, and hence on n, except for values of n 2 
near the peak of the effectiv e po tentia l. T hese consider- 
ations, combined with Eqs. (2.7) and (|2.8|), give 



1 2/(l+3 7i )-2/(l+3 7/ ) 



for inflationary cosmologies, and 



,l+2/(l+3 7< )+2/(l+3 7/ ) 



(2.9) 



(2-10) 



for bounce cosmologies. 

A special case is when 7/ » 1/3 (the equation of state 
for a relativistic or radiation-dominated fluid) , for which 
the effective potential a" /a « 0. This is depicted in 
Fig. ||. If the transition to radiation dominance is rapid, 
then r]f, and hence a/ and a'j, are independent of n for 
a wide range of n. However, waves that are well outside 
the Hubble radius at the momen t of transition will have 
a'f/af 3> n, as is clear from Eq. ( |2.5| ) and Fig. ||. So the 
appropriate terms in Eqs. (2.7) and (2.8) have behaviour 
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This is the same dependence as one woul d get by naively 
plugging 7^ = 1/3 into Eqs. ( pTgj ) and ( |2.10 ), so these 
equations are correct even as 7/ — > 1/3. (This argument 
also applies for waves that hit the effective potential dur- 
ing radiation-dominated collapse.) 
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FIG. 2. Schematic of the evolution of the effective po- 
tential and Hubble scale for inflation (77 < 771) followed by a 
radiation-dominated equation of state (77 > 771). Note that 
rjf = 771 independent of n for a broad range of frequencies. 

A similar situation occurs for waves that remain out- 
side the Hubble radius throughout an intermediate pe- 
riod of 7 « 1/3, when formally they have emerged from 
the effective potential and are oscillatory. This is de- 
picted in Fig. H The effect of an intermediate period 
A77 of zero effective potential is a term ~ sin(rtA77)/n in 
the final amplitude. If the waves remain well outside the 
Hubble radius throughout this era, then 77A77 <C 1, and 



this t erm is a c onsta nt with respect to n. Once again, 
Eqs. (HJ) and ( |2~To| ) remain unchanged. Both this and 
the preceding cases confirm one's physical intuition, that 
the final spectrum cannot resolve the details of sudden 
phase transitions that occur when the waves are much 
larger than the Hubble radius. 




FIG. 3. Schematic of the evolution of the effective po- 
tential and Hubble scale for inflation (77 < 771), followed by 
radiation dominance (771 < 77 < 772), followed by generic decel- 
erated expansion (77 > 772). Waves emerging from the poten- 
tial barrier after 772 will have been outside the Hubble radius 
throughout the period of radiation dominance. 



Eqs. (2.S), ( 2.1C ) also break down when jij w —1/3, 
but this condition will never arise. When 7 w —1/3, one 
has a ~ e K,? , so a" /a ps constant, and waves will neither 
hit nor leave the effective potential. 



B. The initial and final spectra 

Most analyses of primordial gravitational waves as- 
sume that the initial state of the metric perturbations 
(during the Planck era or during pre-Big-Bang collapse) 
is a quantum-mechanical vacuum. While this is a reason- 
able assumption, one does not know the initial state in 
advance of observation. In fact, primordial gravitational 
waves are a potential means of observing the initial state 
of the cosmos. It is instructive, therefore, to leave the 
initial graviton spectrum as a free parameter. 

I parameterize the initial spectrum using the mean 
quantum occupation number of graviton modes as a func- 
tion of mode frequency, N(n). In the semiclassical ap- 
proach, the creation of gravitons can be treated as the 
classical amplification of the vacuum energy hu/2 in each 
mode, so I normalize N(n) to be 1/2 for a mode in the 
ground state (no real gravitons), 3/2 for the first excited 
state (one real graviton), and so on. This approach is 
valid provided the final state is classical (\f3 n \ 3> 1). 

Cosmological perturbation spectra are normally ex- 
pressed as a normalized energy density per logarithmic 
frequency interval, £1 = (dE / dVdlnoj) / p c (where p c is 
the energy density required to make the universe spa- 
tially flat). There are 2oj 2 dujdV modes in a frequency 
interval dio in a volume dV, each with energy hujN(n) 
(where n = aco), so the initial energy spectrum is: 
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Note that ai„itiai is the (constant) scale factor of the uni- 
verse at the time that one is evaluating the initial spec- 
trum; it is not the same as a,, which depends on n. 

In considering the logarithmic slope of the spectrum I 
ignore overall amplitude factors, which are highly model- 
specific. The energy spectrum at the present day is 
quadratic in \(3 n \, so its logarithmic slope is: 



dlnCl _ dlnSljnitiai ^^d\n\(3 n \ 
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where the evaluation is at n 
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^/flp rcscnt- Combining 
this with Eqs. (2.12), (2.9), and (2.10) gives the pricipal 
result of this paper: 
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d In w d In n \ 7$ + 1 /3 , 
for an inflationary cosmology, and: 
dlnn d\nN _ ( 2 7i 



d In w dlnn \ 7, + 1 / 3 
for a bounce-type cosmology. 
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III. COMPARISON WITH SPECIAL CASES 

The formulae given above are consistent with results 
established previously in the literature. I first show this 
for the two "standard" cases of deSitter inflation and 
of cosmological rebound from a state of weakly coupled 
strings. I then compare with recent papers that consider 
more generic pre- and post-inflationary evolutions. 

A. Standard inflation 

The "standard" model for an inflationary cosmology 
consists of a period of deSitter expansion (7 = —1) due to 
the universe being in a false vacuum with nonzero energy 
density. In the post-inflationary era, the universe quickly 
reheats to a radiation-dominated equation of state (7 = 
1 /3), followed by cooling to a matter-dominated equation 
of state (7 = 0). If we assume the initial graviton state 
from the Planck era was a vacuum, then for waves that 
emerged into the radiation-dominated universe we have 
N = 1/2, 7l = -1, 7/ = 1/3, and 



d\nn 
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(3.1) 



This is, of course, the flat spectrum predicted by Harrison 
and Zel'dovich ' 

The last two decades of the spectrum (10~ 16 - 
10 -18 Hz) consist of waves that emerged more recently, 
when the universe was matter-dominated (7/ = 0). In 
this band the spectrum is tilted towards low frequencies, 
with a spectral index of —2. 



B. Standard string- motivated rebound cosmologies 

It has been shown |||| that string theories can lead to 
a gravitational-wave background that is increasing with 
frequency. These models are consistent with the cur- 
rent derivation provided one recognizes that the current 
formalism applies to the Einstein frame (where the met- 
ric represents intervals measured by physical rulers and 
clocks), rather than the string frame (where the met- 
ric geodesies describe the trajectories of weakly coupled 
strings). In the Einstein frame, a typical evolution con- 
sists of an epoch of decreasing lengthscale a and weakly 
coupled strings, followed by a period of strongly coupled 
strings that reverse the collapse, followed by relaxation of 
the string dilaton and an expanding radiation-dominated 
universe, eventually cooling to matter dominance. As 
usual, a vacuum initial state {N = 1/2) is normally as- 
sumed. Most modes hit the effective potential during 
the epoch of weakly interacting strings, when the kinetic 
term of the dilaton dominates the energy background, 
giving a nearly maximally stiff equation of state (7, w 1). 
For waves that emerge dur ing the radiation-dominated 
phase (7/ w 1/3), Eq. ( 2.15 ) gives: 



din ft 
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Eq. (3.3) of H and Eq. (5.7) of |j) also give a spectral in- 
dex of 3 for low frequencies. These papers also predict a 
logarithmic cutoff at high frequencies, ~ uj 3 In 2 {uj s j u), 
where w s is the maximum frequency of waves that en- 
countered the effective potential during the epoch of max- 
imal stiffness (typically of order 10 11 Hz, depending on 
the details of the strongly coupled string epoch). This 
effect occurs only for 7, exactly equal to 1, which the 
current analysis doesn't consider, and only affects the 
spectral index within a decade of tu s in any case. 

The low- frequency tail (ui < 10~ 16 Hz) of waves that 
emerge when 7/ = will also be tilted towards high 
frequencies, but with a spectral index of 1. 



C. Nonstandard equations of state 

A recent paper [0 has explored the effects of an inter- 
mediate phase between the inflationary and radiation- 
dominated eras, in which the dominant cosmological 
fluid had an equation of state stiffer than radiation 
(1/3 < 7/ < 1). One consequence is that the background 
gravitational-wave spectrum would be tilted towards 
high frequencies for those wavelengths that emerged dur- 
ing this era. Assuming an in itial vacuum (N = 1/2) and 
deSitter inflation, Eq. (2.14) gives us: 



1/3 



dlnu 7/ + 1/3 
Eq. (3.32) of M gives the same dependence on 7. 



(3.3) 
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Eq. (3.31) of J5J extends the analysis to the case of 
a maximally- stiff equation of state (7/ = 1), which this 
paper does not consider. This gives a high-frequency log- 
arithmic cutoff Q ~ lu In 2 (uji/lu), where w\ is the high- 
est frequency that encountered the amplifying potential. 
The cutoff does not affect the spectral index at frequen- 
cies more than a decade or so below cji, which is likely 
in the MHz to GHz range (well above the pass-bands 
of proposed gravitational- wave detectors) . This is anal- 
ogous to the situation mentioned earlier, when 7; -> 1 
during collapse. 

Gasperini Q has performed an even more general anal- 
ysis of the primordial gravitational wave spectrum, con- 
sidering generic evolution of 0,(77) both in the early and 
late cosmological epochs. As usual, the initial state of th e 
waves was taken to be a vacuum (N = 1/2). Eqs. (2.14), 
( [2.15 ) can then be written as: 



or 
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for monotonic or bounce cosmologies, respectively. By 
comparison, Eq. (4.26) of gives the spectral index as 
4 — 2^1 — 2^+1, where v = \a — 1/2| and a is the exponent 
of the power law a ~ rj a during the epoch when the waves 
first strike the effective potential (for fj+i) and when 
they leave the effective pot enti al (fo r W ). Gasperini's 



results are identical to Eqs. (5A) and (3.5), provided one 
recognizes two things. 

First, a phase of slowly varying a corresponds to a 
phase of slowly varying 7, with a = 2/(1 + 37). In- 
flationary cosmologies have ji < —1/3, hence ai+i < 0. 
Bounce cosmologies have — 1/3 < 7, < 1 and ttj+i > 1/2. 
Either cosmology has —1/3 < 7/ < 1 and ai > 1/2. (The 
condition < a < 1/2 corresponds to the physically un- 
realistic case of 7 > 1.) 

Second, Gasperini assumes in Eq. (4.26) of H that all 
waves leave the effective potential at the same cosmolog- 
ical epoch, so that v\ is a constant across all frequencies 
(normally v\ = 1/2 for radiation dominance). In fact, 
the spectral index (but not necessarily the normaliza- 
tion) given in Eq. (4.26) of W is valid for v\ varying 
across frequencies. 



IV. CONCLUSIONS 



Eqs. ( |2.14| ) and ( 2.15 ) show how much can be deter- 
mined about early cosmology from an observation of the 
spectral index of primordial gravitational waves. Al- 
though these formulae are quite simple, they still involve 
three independent parameters: 7$, 7/, and N(n). To 
make definite statements about any one of them, one 
must make assumptions about the others. The usual 



procedure is to assume that we know N(n) and 7/, leav- 
ing 7^ to be deduced; however, it is quite possible that 
future advances in theory will fix both ji and 7/, allow- 
ing gravitational-wave observations to probe directly the 
Planck-scale structure of the universe. 

The formulae are also useful in showing which cosmo- 
logical theories produce equivalent gravitational signa- 
tures. For instance, a monotonically inflating cosmology 
with 7^ < — 1/3 will yield the same spectral index as a 
bounce-type cosmology with 7i 2 = (7^ + 1)/(37; 1 — 1) 
in the range (—1/3, 1/3). Since these are physically valid 
ranges for 74 in each case, it follows that a monotonic cos- 
mology cannot be distinguished from a bounce cosmology 
given only the gravitational-wave spectral index. 

The purpose of this paper was to present the sim- 
plest possible expressions relating a generic cosmologi- 
cal gravitational- wave signature to the physical processes 
that produced it, with a minimum of assumptions about 
the underlying cosmological model. I have deliberately 
chosen an approach of "maximum ignorance" , allowing, 
as much as possible, for the observations to dictate the 
cosmological parameters. Formulae such as these may 
prove valuable in an era when gravitational-wave obser- 
vations begin to explore a new realm of physics, possibly 
revealing things beyond our most informed predictions. 
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